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Введение. Проблемы динамики процесса резания, рассматри- 
ваемые в единстве устойчивости равновесия в подвижной 
системе координат, определяемой траекториями исполни- 
тельных элементов станка, а также свойств притягивающих 
множеств, формируемых в окрестности точки равновесия, 
являлись предметом многочисленных исследований. В отли- 
чие от этих работ, в данной статье обсуждается проблема 
формирования притягивающих множеств, формируемых в 
динамической системе резания, в которой траектории движе- 
ния исполнительных элементов станка являются возмущен- 
ными периодическими колебательными смещениями шпин- 
дельной группы. Эти возмущения во многом определяются 
биениями шпиндельной группы. Они вызывают периодиче- 
ские колебательные смещения инструмента относительно 
обрабатываемой детали и в совокупности с динамическими 
свойствами системы приводят к различным динамическим 
эффектам. Во-первых, за счет параметрического самовозбуж- 
дения возможна потеря устойчивости установившихся стаци- 
онарных траекторий периодических деформационных смеще- 
ний. Во-вторых, в зависимости от частоты и амплитуды про- 
странственных биений шпиндельной группы, возможны фор- 
мирования притягивающих множеств типа предельных цик- 
лов и инвариантных торов и их вырождение. Наконец, при 
вариациях частоты имеют место бифуркации удвоения перио- 
да колебаний и, как следствие, образование хаотических ат- 
тракторов. Причем их формирование является чрезвычайно 
чувствительным, как к параметрам взаимодействующих через 
процесс резания подсистем со стороны инструмента и обраба- 
тываемой детали, так к параметрам динамической связи фор- 
мируемой процессом резания. 

Материалы и методы. Свойства притягивающих множеств и 
законы их формирования изучаются на основе цифрового 
моделирования. Для этого приводится математическая модель 
возмущенной системы, а также результаты изучения форми- 
рования притягивающих множеств в динамической системе 
резания на примере процесса точения. 

Результаты исследования. Приводятся примеры влияния 
биений на притягивающие множества системы, в том числе 
рассматриваются бифуркации системы в зависимости от па- 
раметров возмущений. 
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Постановка проблемы. Синергетический подход к анализу и синтезу управляемого процесса обработки на 
станках полагает рассмотрение двух взаимодействующих через зону резания подсистем со стороны инструмента и 
обрабатываемой детали [1-5]. В отличие от традиционных методов, управление в этом случае формируется на основе 
согласования внешнего управления (программы ЧПУ) с динамическими свойствами процесса обработки [2—5]. Необ- 
ходимо учитывать не только устойчивость траекторий [6-8, 12, 13, 15], но и формируемые притягивающие множества 
(предельные циклы, инвариантные торы, хаотические аттракторы и пр. [9-—11, 14]). При построении математической 
модели динамической связи между подсистемами раскрывается зависимость сил резания от траекторий исполнитель- 
ных элементов станка и упругих деформаций вершины инструмента относительно детали несущей системы станка 
[6-17]. В этом случае учитываются экспериментально наблюдаемые эффекты зависимости сил резания от скорости, 
их запаздывание по отношению к изменениям смещения инструмента относительно детали, наличие следа на детали 
от траектории на предыдущем обороте и пр. [16, 17]. Ранее авторами были рассмотрены проблемы влияния флуктуа- 
ций движения исполнительных элементов станка на траектории формообразующих движений вершины инструмента 
относительно детали [18—20] и на их устойчивость [21]. Настоящая статья является развитием основных положений 
ранее опубликованных работ [18—21]. Анализируются закономерности изменения свойств притягивающих множеств в 
динамической системе резания, которые вызваны флуктуациями параметров динамической связи, обусловленные би- 
ениями шпиндельной группы. Например, если без биений шпинделя в системе формируется притягивающее множе- 
ство типа предельного цикла, то периодические смещения инструмента относительно детали, вызванные биениями 
шпиндельной группы, приводят к взаимодействию изменения параметров системы с предельным циклом. При этом в 
зависимости от траектории деформационных смещений, вызванных биениями, могут наблюдаться различные бифур- 
кации притягивающих множеств (предельных циклов, двумерных торов, хаотических аттракторов). Необходимо учи- 
тывать, что именно тип и параметры притягивающих множеств непосредственно влияют на топологию формируемой 
резанием поверхности и ограничивают возможности синергетического управления процессами обработки на металло- 
режущих станках. 

Математическая модель. Флуктуации параметров, рассматриваемые в статье, обусловлены биениями шпин- 
дельной группы станка, зависящие от биений шпинделя, геометрических погрешностей изготовления и установки за- 
жимного приспособления, несовпадения оси вращения шпинделя и оси задней бабки и др. (рис. 1). Они являются пе- 


риодическими функциям времени ДАХ, (Г) = АХ, (Е—АТ), К=1,2,3...,1=12,3. Смещениям АХ, (Г) соответствуют ско- 
рости АТ (1) = 4АХ, (1 / 4Ё. Притягивающие множества, формируемые при резании в возмущенной системе, образуют- 


ся на основе взаимодействия флуктуаций динамической системы со свойствами невозмущенной системы. Для рас- 


крытия закономерностей влияния АХ, (1), АТ (1) на притягивающие множества ограничимся случаем обработки абсо- 


лютно жесткой детали. Тогда упругая система определяется только подсистемой инструмента, в которой деформации 
отсчитываются от баз станка, привязанных к его несущей системе (рис. 1). Для описания наиболее важных свойств 
подсистемы инструмента она схематизируется в виде элемента с обобщенной массой, подвешенного в пространстве 
на упруго-диссипативных подвесках [9-11]. Биения заготовки приводятся к вариациям точки контакта вершины ин- 


струмента без учета упругих деформаций детали. 
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АХ (1) ={АХ 1, АХ 2. 
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Рис. 1. Ориентация осей деформационных смещений, возмущений и сил, действующих на режущий инструмент 
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Динамику такой системы можно описать уравнением, обоснование которого дано в работах [23-25]: 


2 
Хх Хх 
И а (1) 
(И (ИА 
ле Х={Х,Х,,Х . — вектор упругих деформационных смещений вершины инструмента в [мм]; 
т=[ тк т, =т, при:5=К^, т,, =0, при:5 К, 5,К=1,2,3 в кгс? /мм , й=|[ И, , |5. К =1,2,3 в  кКгс/мм, 


с = Е. =1,2,3 в кг / мм — симметричные и положительно определенные матрицы инерционных, скоростных и 


упругих коэффициентов подсистемы инструмента; Р = {[Н +Ео о |, в, +Ро — силы резания в координатах 


ыы 1 а к 
СОСТОЯНИЯ И внешних воздеиствиях. Здесь силы Ё | 5 р _ ‚ деиствующие на задние грани, зависят от скоростеи их 


сближения с обработанной частью детали и длины контакта (Г и [3 на рис. 1) рассматриваемой грани с поверхностью 


заготовки. То есть они представимы в виде 


о. =р.{ | [Из —аХ, /аЦ- АХ, (}ехр(АИ (1) - си аХ, / 40); С) 
-Т 
ЕЗ, = р. [1 +АХ (1) Х, (В ]ехр, , — 0-Х, / &+ АУ, (0], 


где р, р. — удельные значения сил, приведенные к длине контакта задних граней инструмента с поверхностью об- 
1 Е. 
работанной части детали в [кг / мм]; 9, @&, — коэффициенты крутизны нарастания сил по мере сближения поверхно- 


стей инструмента и детали. 
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ТАК /@+Е = рож, + рехр[-окИ, › + АУ, (ах, Чи + 
АХ О- ХО, АХ, 1+ АР, (48 
ТГ 


Т.Е, /@+Е, = рох› А+ нехр[-оКИ, › + АТ, (0) -аХ, / 41] [1ь” + 
| 3 
+АХ, (1) - Х, (#]] [| {7,. -аХ, Гаваё (3) 
1-Т 
ТЕ, /@+Е. = роз + нехр[-о(И, ‚ + АУ, (1) -аХ, Г} 1” + 


+АХ (0-Х, АХ, аа, 


Здесь Л; при! =1,2,3 — силы, действующие на переднюю грань инструмента, с учетом запаздывания, определяемого 
постоянными времени Г,1=1,2,3; р, — давление стружки на переднюю поверхность инструмента в области малых 
скоростей резания; и — коэффициент, определяющий соотношение сил резания в областях малых и больших скоро- 
стей; о — коэффициент, определяющий убывание сил по мере увеличения скорости резания; И, , И. — 
определяемые технологическими режимами значения скоростей резания и продольной подачи; Г — время оборота 


детали. Так как 1 = 1, = 1., то в переходных процессах учитывается изменение ориентации сил в пространстве. В 


х Т 
установившемся состоянии без возмущений ориентация сил неизменна и определяется коэффициентами {).,х., Хх.) 


Модель (1), (2), (3), как и ранее [13], учитывает: уменьшение по экспоненциальному закону сил при увеличении ско- 
рости, запаздывание вариаций сил по отношению к деформационным смещениям, закономерности формирования но- 
вой площади срезаемого слоя. 

Эта модель отличается тем, что учитывает изменение ориентации сил в нестационарном режиме, формирова- 
ние дополнительных сил за счет сближения задних граней инструмента с обрабатываемой деталью в. . о Эти 
силы зависят не только от свойств сближения задних граней с поверхностями заготовки, но и от изменяющейся в ди- 
намике длины контакта инструмента с заготовкой. Это существенное обобщение и уточнение известных положений о 
неоднозначности изменения сил, действующих на инструмент, при его сближении с заготовкой и отталкивании [22]. 
Эти свойства характеризуют нелинейную диссипацию процесса резания. В ней приняты также во внимание возмуще- 
ния АХ, (1), АЙ (1) ‚ как заданные периодические функции АХ (1) = АХ (1-АТ), АГ( =АТ(Е-АТ,). 

Бифуркации притягивающих множеств. Изучение бифуркаций притягивающих множеств выполнено на 


основе цифрового моделирования системы (1-3) в зависимости от параметров биений. Биения представляются в виде 
АХ (В = АХ, {эт ОБзт(ОЕ+л/2)}', то есть рассматриваются в плоскости, нормальной к оси вращения детали. Траек- 
тории биений представляют окружности. Преобразование притягивающих множеств выполнено по следующей схеме. 
Рассматривался процесс врезания инструмента и образования стационарных состояний без биений шпинделя. Затем 
изучались переходные процессы и установление нового стационарного состояния при введении в систему биений. 
Диаграмма бифуркаций притягивающих множеств выполнена по сетке, при этом было вычислено более 10 000 соче- 
таний исходных состояний системы, параметров биений и динамической связи, формируемой процессом резания. 


Проанализирована система, параметры которой приведены в табл. 1, а модели сил в координатах состояния в табл. 2. 
Обобщенная масса — т= 0,035 кг / с”. Для удобства принято Т. = Г, =0,5Т, = Т,. Для этого случая ранее проанализи- 
рованы механизмы потери устойчивости траекторий [21]. Примеры вычислены для режимов: продольная подача 
5р = (0,05—0,3) мм, глубина &, =(0,5—2,5)мм ‚ скорость Г, „ = (500-2000) мм/ с. Связь параметров модели с тех- 
нологическими режимами изложена в более ранних публикациях [2, 12-14, 16, 17]. Если свойства системы без биений 
заданы, то их преобразования, вызванные биениями, зависят от амплитуды АХ. и частоты &). Здесь наибольший ин- 


терес представляет ситуация, когда без биений в системе формируется орбитально асимптотически устойчивый пре- 


дельный цикл, то есть автоколебания. 
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Таблица 1 
ТаЫе 1 
Параметры динамической системы 


Дупатгс 5уЯяет рагатеету 


пт, й, ›, й, › — |, Й, 3 —- й т, Й, 3 — й 2, 
кг.с/ мм кг.с/ мм кг.с/ мм кг.с/ мм кг.с/ мм 
2 2 


ОО ОЕ ПОСЕВ 





Таблица 2 
Табе 2 
Параметры модели сил в координатах состояния 


Рагатеегу оГюгсе тоае[ т ие соот@тае$ 





(50—900) (0,01—0,1) (0,2-0,5) (0,2-0,5) (0,01-0,1) 


Частотный диапазон биений можно разделить на три поддиапазона. Низкочастотный диапазон (в пределах 
0,1-10,0 Гц), в котором вдоль периодической траектории параметры системы можно считать медленно изменяющи- 
мися, а систему анализировать как замороженную. Среднечастотный — в пределах 10,0—400,0 Гц, в котором происхо- 
дит сложное взаимодействие автоколебаний с колебаниями, обусловленными параметрическими возбуждениями. Вы- 
сокочастотный — выше 400,0 Гц, в котором частоты биения находятся за пределами полосы пропускания колебатель- 
ных контуров. Однако и в этом случае биения влияют на динамическое качество за счет эффектов вибрационного 
сглаживания нелинейных связей, формируемых резанием. 

Вначале рассмотрим низкочастотный поддиапазон. На рис. 2 приведен пример динамической перестройки си- 
стемы при ее возмущении гармоническими биениями шпинделя. На рис. 2а звездочкой отмечена точка, для которой 
приведены примеры траекторий. 
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Рис. 2. Пример динамической перестройки системы при движении вдольвозмущенной траектории с частотой 0,3 Гц: а — область 
устойчивости невозмущенной системы в плоскости «ру — 1ь »; 6 — установление стационарного периодического режима; 


с — увеличенный фрагмент динамической реорганизации деформационных смещений на участке «Е-Е» 


Ето. 2. Ехатре оГ5уяет Аупат:с аиятет ийеп тоуте аопе ретигбей гадесогу иийй [гедиепсу ог 0.3 Н2: а- 5аБИйу 
гетоп орипрепитЬе4 зуяет т “ру — Ть ”рапе; Ь — зейтз о} байопагу репо с тоде; с — тсгеазе4 Ггавтет оГ4упат!с 


геатгапеетет! оГаеютгтайоп @5р1асетеп т “Е-Е” 5еспоп 


На траекториях выделены участки: «0-А» — переходной процесс врезания инструмента в деталь; «А-В» — 
стационарный режим без биений; «В-С» — переходной процесс установления стационарного периодического режима 
с биениями. В выделенной точке равновесие является неустойчивым, в системе формируется притягивающее множе- 
ство типа предельного цикла (на иллюстрации выделено в окружности). Частота автоколебаний в этом случае без би- 
ений равна 175 Гц. При возмущении системы с частотой биений 0,3 Гц вдоль траектории изменяются параметры си- 
стемы, так как меняется текущее значение припуска и скорости деформационных смещений в плоскости. Однако это 
медленные изменения, которые вызывают изменения уравнений в вариациях относительно рассматриваемой точки. 
Наблюдается эволюционная перестройка притягивающих множеств (на участке «С-О») и свойство устойчивости (на 
участке «О-Е» траектория является асимптотически устойчивой). 

По мере увеличения частоты биений в пределах от 10,0 Гц до 400,0 Гц в системе, в зависимости от частоты и 
амплитуды, наблюдаются следующие эффекты нелинейных взаимодействий (рис. 3). При малой амплитуде биений в 
установившемся состоянии наблюдается формирование двумерного инвариантного тора. При дальнейшем увеличении 
амплитуды имеет место преобразование тора в предельный цикл на частоте возмущений. При этом имеет место эф- 
фект смещения упругих деформационных смещений инструмента, усредненных по периоду колебаний, зависящего от 
амплитуды и частоты биений шпинделя (рис. 4). Он связан с тем, что нелинейные связи процесса резания не обладают 
свойством центральной симметрии относительно биений. Это, прежде всего, связи, формируемые контактом задних 
граней инструмента с обработанной частью детали. Таким образом, текущее значение диаметра обрабатываемой дета- 
ли зависит не только от сил резания и параметров упругости взаимодействующих подсистем, но и от параметров 
установившихся периодических движений. Этот эффект особенно заметен в том случае, если частота биений близка к 
частоте автоколебаний системы без биений (на участке «А-В» на рис. 4). В общем случае при варьировании частоты и 
амплитуды биений шпинделя наблюдаются эффекты формирования двумерного инвариантного тора, его вырождения 
в предельный цикл и образования хаотических аттракторов, который формируется после каскада бифуркаций удвое- 
ния периода колебаний (рис. 3). Однако формирование хаотических аттракторов наблюдается не во всех частотных 
диапазонах. 
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Рис. 3. Диаграмма бифуркаций притягивающих множеств: 
| — вынужденные одночастотные колебания (предельный цикл); 
2 — двумерный инвариантный тор; 
3 — установление хаотического аттрактора через каскад удвоения периода; 
4 — система, неустойчивая в целом 


Ето. 3. ДЛаэтат огайгасиие 5ез Буигсаноп5: 
1 - югсе4 зта[е-гедиепсу озсШайоп5 (ти сус[е); 
2 - 1ло-@тепяопа![ туатапЕ югиб; 
3 - вет оГсйаойс айгасюг йгоией репо4-адоиб Пиз сазсаае; 
4 - вбуяет ипяаЫе а5 а уйо[е 


Во всех случаях при фиксированной частоте существует предельное значение амплитуды биений, начиная с 


которой система становится неустойчивой в целом (область 5 на рис. 3). Причем, по мере увеличения частоты биений 


шпинделя, критические значения амплитуды периодически изменяются. Кроме этого, имеется тенденция уменьшения 


критической амплитуды по мере увеличения частоты. 
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Рис. 4. Пример переходных процессов деформационных смещений инструмента при частоте радиальных биений 


шпинделя равной 180 Гц: а — временные траектории; Б — фазовая траектория для Х’.На иллюстрации 
| 


(АХ! ,, АХ, ,, АХ, ,) — динамическое смещение точки равновесия подсистемы инструмента 
Ето. 4. Ехатре отапяеп! ргосез5ез оГоо[ аеогтаноп @5р[асетепи а! гедиепсу орзртаЁе га Беайиз едиа![ ю 
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Особенно наглядно эти эффекты можно проследить, если рассмотреть фазовые траектории деформационных 
смещений. Приведем некоторые примеры (рис. 5). Они соответствуют точкам, выделенным на рис. 3. Как видно, по 
мере увеличения амплитуды биений имеют место следующие преобразования притягивающих множеств. Начальному 
состоянию системы без биений шпинделя соответствуют автоколебания. При малых амплитудах биений шпинделя 
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наблюдается простое наложение автоколебаний на колебания, определяемые изменениями параметров, связанными с 
биениями. При дальнейшем увеличении амплитуды имеет место затягивание траектории к частоте, соответствующей 
частоте биений. Таким образом, инвариантный двумерный тор преобразуется в предельный цикл. Дальнейшее увели- 
чение амплитуды вызывает каскад бифуркаций удвоения периода и формирование хаотического аттрактора. При 
дальнейшем увеличении амплитуды наблюдается потеря устойчивости системы в целом. 


АХ /армм/с АХ 1 / 4, мм/с АХ / @,мм/с 






13 13,5 14 14,5 Х,х10?, мм т 20 
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Рис. 5. Преобразование фазовых траекторий в плоскости ( А" —ЯХ\ / Е) для установившегося состояния 
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При увеличении частоты от 400 Гц и выше наблюдается практически полное наложение вынужденных коле- 
баний на формируемые без биений автоколебания. При этом переходной процесс при установлении стационарного 
режима, соответствующего двумерному тору, становится малым, и он уменьшается по мере увеличения частоты бие- 
ний шпинделя. Такое свойство наложения биений на автоколебания объясняется тем, что в рассматриваемом частот- 
ном диапазоне частоты биений лежат за пределами полосы пропускания упругой системы инструмента. Колебатель- 
ные же смещения, обусловленные биениями, непосредственно прикладываются к режущему инструменту. 

Анализ результатов. Свойства динамической системы резания принципиально зависят от динамической свя- 
зи, формируемой процессом резания, характеристик взаимодействующих подсистем и вибрационных возмущений, 
определяемых биениями шпиндельной группы станка. Традиционные представления о преобразовании биений шпин- 
дельной группы непосредственно в траектории формообразующих движений инструмента относительно детали и, 
следовательно, в характеристики геометрического качества формируемой резанием детали, далеко не всегда соответ- 
ствуют действительности. Как показал анализ, установившиеся колебания деформационных смещений вершины ин- 
струмента относительно детали образуют различные притягивающие множества: предельные циклы, инвариантные 
торы и хаотические аттракторы. Причем, конфигурация, например, предельных циклов может быть достаточно слож- 
ной, так как вдоль траектории периодических движений изменяются свойства динамической связи. Более того, по ме- 
ре увеличения частоты биений шпиндельной группы, вызывающие периодические изменения параметров системы, 
имеют место различные параметрические явления, в том числе параметрическое самовозбуждение системы. Кроме 
этого обнаружены эффекты динамического смещения точки равновесия, низкочастотной модуляции этого смещения. 
При определенных параметрах динамической связи и свойствах биений формируется хаотическая динамика. Как пра- 
вило, она формируется после каскада бифуркаций удвоения периода колебаний. Все эти факторы, изменяя траектории 
формообразующих движений вершины инструмента, не только усложняют прогнозирование параметров точности и 
качества изготовления деталей, но и приводят к непредсказуемым эволюционным их изменениям. 
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Заключение. Периодические флуктуации параметров динамической характеристики процесса резания во 
многом определяются пространственными периодическими биениями шпиндельной группы с учетом точности уста- 
новки заготовки. Флуктуации не только изменяют траектории формообразующих движений, но и влияют на динами- 
ческое качество процесса резания, в том числе устойчивость траекторий. Приведенный материал, созданные матема- 
тические модели и методики изучения влияния возмущений на устойчивость не только дополняют знания о динами- 
ческой системе резания, но и позволяют оценивать предельно достижимую точность обработки при управлении про- 
цессом резания и создании программы ЧПУ при обработке деталей на конкретном станке. Таким образом, для каждо- 
го станка существуют предельно достижимые показатели качества изготовления деталей, диктуемые текущим состоя- 
нием станка. Приведенный материал позволяет также оценивать достижимые параметры качества обработки на кон- 
кретном станке по параметрам биений шпинделя. 
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